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Abstract 

In this paper we investigate the regularizing behavior of two-phase 
Stefan problem near initial data. The main step in the analysis is to 
establish that in any given scale, the scaled solution is very close to a Lip- 
schitz profile in space-time. We introduce a new decomposition argument 
to generalize the preceding ones ( CJK1 - CJK2 and CK ) on one- phase 
free boundary problems. 

1 Introduction 

Consider u {x) : B R (0) — > 1R with R » 1 and u > —1, \{u = 0}| = and 
uo(x) = —1 on dB R (0). (See Figure 1.) The two-phase Stefan problem can be 
formally written as 

u t - Au = in {u > 0} U {u < 0} 

\Du + \ - \Du~\ on d{u>0} 



(ST2) 



u(-,0) = u 

u = -l on dB R (0) 

where u + and u~ respectively denotes the positive and negative parts of u, i.e, 

u + := max(u,0) and u~ := — min(w, 0). 

The classical Stefan problem describes the phase transition between solid/liquid 
or liquid/liquid interface (see [M] and also |OPRj .) In our setting, we consider a 
bounded domain O C Br(0) and the initial data u (x) such that {u > 0} = SIq 
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Figure 1: Initial setting of the problem 



and {uq < 0} = M n — Qq. To avoid complicity at the infinity, we consider the 
problem in the domain Q = -Br(O) x [0,oo), with Dirichlet condition 

u = f(x,t) < on dB R (0), 

where f(x,t) is smooth. In (ST2) we have set / = — 1 for simplicity. Since 
our initial data will be only locally Holder continuous, we employ the notion of 
viscosity solutions to discuss the evolution of the problem. Viscosity solutions 
for (ST2) is originally introduced by }ACSlj (also see }CSj ). As for existence 
and uniqueness of viscosity solutions, we refer to jKPj . 

Note that the second condition of (ST) states that the normal velocity V x .t 
at each free boundary point (x, t) £ d{u > 0} is given by 

V x , t = \Du+\ - \Du-\(x,t) = (Du+(x,t) - Du-{x,t)) ■ v XtU 

where v x t denotes the spatial unit normal vector of d{u > 0} at (x, t), pointing 
inward with respect to the positive phase {u > 0}. 

In this paper we investigate the regularizing behavior of the free boundary 
d{u > 0}. Our main result states that when Tq has no sharp corner, then 
the free boundary immediately regularizes after t = 0, and stays regular for a 
small amount of time. Note that, in general, after some time the free boundary 
may move away from its initial profile and develop singularities by topological 
changes, such as merging of two boundary parts. Whether this happens with 
star-shaped initial data is an open question (see Remark |3. 11 ) 

The well-known results of }ACSlj - fACS2j states that if a solution as well 
as its free boundary of (ST2) stays close to a locally Lipschitz profile in a 
unit space-time neighborhood, then the solution is indeed smooth in a smaller 
neighborhood. Hence the main step in our analysis is to prove that the free 
boundary d{u > 0} stays close to a locally Lipschitz profile over a unit time 
interval. Indeed proving this step has been the main challenge in the previous 
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work of the authors and Jerison ( jCJKlj . |CJK2j . jCKj ) on the study of one- 
phase free boundary problems. Once this step is established, using the fact that 
u is a caloric function in almost Lipschitz domain, we will have some control over 
the behavior of re-scaled solutions following the arguments in [CJKlj . Then the 
appropriate modification of iteration arguments taken in |ACS1] - |ACS2] applies 
to derive further regularity results (see section 5). In extension of the ideas 
from one-phase case to our setting, the main obstacle lies in the competition 
between fluxes of positive and negative phase: to overcome this, we introduce a 
new decomposition procedure which we explain below. 

Before discussing our result in detail, let us introduce precise conditions on 
the initial data. 

(I-a) fio aud uq are star-shaped with respect to a ball B ro (0) C f^o- 

Observe that then the Lipschitz constant L of OQq is determined by r and 
do, where da := svop{d(x, B ro (0) : x <G 8Q,q}: i.e., there exist h = h(ro) and 
L = L(^) such that for any xo <G dflo, after rotation of coordinates one can 
represent 

B h (x ) n n Q = {(x', x n ) : x' G M n ~\ x n < f(x)} with Lipf < L. (1.1) 

For simplicity of the presentation we set h = 1. 

For a locally Lipschitz domain such as fto, there exist growth rates < (3 < 
1 < a such that the following holds: Let if be a positive harmonic function in 
^oHi?2(a;), x £ <9r2o, with Dirichlet condition on d£lo^B2{x) 1 and with value 1 
at x — e n . (Here let e„ be the direction of the axis for the Lipschitz graph near 
x.) Then for x — se n €E Oo H B\{x) 

s a < H{x - se n ) < s P . (1.2) 

Now wc precisely describe the range of the Lipschitz constant L of Q . 

(I-b) L < L n for a sufficiently small dimensional constant L n so that 

5/6 < P < a < 7/6. 

The remaining conditions are on the regularity of Uq. 
(I-c) -JV < < N in fi U (B R (0) - O ), 

(I-d) For x G <9O , e„ = x/\x\ and small s > (for < s < 1/10), 

\Du (x±se n )\ > Cs a -\ 
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Note that (I-c) and (I-d) holds for uo which is smooth in its positive and negative 
phases and is harmonic near the initial free boundary: i.e., — A.uq = in the set 
({u > 0} U {u < 0}) n {x : d(x, dn Q ) < 1}. 

For a function u(x,t) : M n x [0,oo) — > M, let us denote 
Q(u) := {u > 0}, il t (u) := {u(-,t) > 0} 

and 

r(u) := d{u > 0}, r t (u) := d{u(-, t) > 0}. 

Since Tq = d{u(-,0) < 0} in our setting, the property is preserved for later 
times, i.e., 

r t (u) = d{u(-,t) > 0} = d{u(-,t) < 0} for alU > 

(see [RB], EZ], and [KP]). 

In the one-phase case the problem can be written as follows: 

' u t - Au = in {u > 0} 



(ST1) 



\Du\ on d{u > 0} 



\Du\ 
I u(-,0) = u+. 
In |CKj the following has been proved for (ST1). 

Theorem 1.1. ([CKJ Theorem 0.1.) Suppose u is a solution of (ST1) in 
-82(0) x [0,1], € Tq(u), with the initial data uq > satisfying (I-b), (I-c) 
and (I-d) in -62(0). Suppose uq{— e„) = 1. If sup B2 (Q) X r 0)1 i u < Mq, then there 
exists a small s > depending on Nq, Mq and n such that the free boundary 
F t (u) becomes smooth and averages out in B s (0). More precisely, 

(a) The free boundary Tt(u) is C 1 and is a Lipschitz graph with respect to e n 
with Lipschitz constant L' < L n in B s (0) C IR n . 

(b) The spatial normal ofT t (u) is continuous in space and time, in B s (0). 

(c) If xe T (u) n B 8 (0) and x + de n E T t (u) n B s (0), then 

x u(x-de n ,0) u{x-de n ,0) 

C < \Du{x + de n ,t) \= V x+ de n ,t < O 

a d 

where C depends on n and Mq ■ Hence 

d u(x-de n ,0) 

- - \Du(x + de n ,t)\ . 

t d 
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Theorem 11.11 states that the free boundary regularizes in space, in a scale 
proportional to the distance it has traveled. Note that the regularity results 
hold up to the initial time and all the regularity assumptions are imposed only 
on the initial data. 

Our aim in this paper is to extend the above theorem to the two-phase 
case. Here the intuition is rather straightforward, based on the previous results. 
There arc two cases: 

(a) One of the phases has much bigger flux than the other: in this case one- 
phase like phenomena (rcgularization by the dominant phase proportional 
to the distance the free oundary traveled) is expected. 

(b) Both phases are in balance: in this case one expects regularization due to 
competition between two phases, resulting in Lipschitz-likc behavior over 
time. 

The difficulty in making above heuristics rigorous lies in introducing a proper 
"sorting" procedure which divides the cases (a) and (b) in a given scale. To 
enable such procedure, it is essential to show Harnack-type inequalities for so- 
lutions of (ST2) in both cases, ensuring that the behavior of solutions can be 
localized in a proper time-space scale. 

To state the main result, we introduce one more notation. 

• For xq £ To = ro(it) and e n := xo/l^ol, define 

D+/ t\ u + (x - de n ,0) u-(x o + de n ,0) 

R + (x ,d) := — — — — , R (x ,d):= — — 

u (xo + de n ,U) u + (xo—de n ,0) 

and 

t{x , a) := 



u + (xq — de n , 0) ' u (xo + de n ,0) 



Theorem 1.2 (Main Theorem). Suppose u is a solution of (ST2) with ini- 
tial data Uq satisfying (Ia)-(Id) with f2o(u) C i?2(0). Then there exists a 
constant do depending on the dimension n and Nq such that the following 
holds. If xq G To(w) and d < do, then T(u) is a Lipschitz graph in space- 
time in the region i?2d(^o) x [*(^0; d)/2, t(xo, d)] with T(u) intersecting with 
Bd(xo) x [t(xo,d)/2,t(xo,d)]. Further, there exists a positive dimensional con- 
stant M such that the following holds. 

(a) IfR+(x ,d) > M, then 

\Du+\(x,t)~ U+{X0 - den > 0) 

. . ,t(xo,d) 
in B2d(Xo) x [ ^ 'H x 0,«jj and 

u + {x Q -de n ,0) 
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on T(u) n (B 2 d{x ) x [ ^ X0 ' d \ t(x ,d)}). 
The parallel statements hold for u~ if R~ (xo,d) > M. 
(b) IfR+{d),R-(d) <M, then 

|2*. |(x,t) - d - d 

in B 2 d(xo) x [t(x ,d)/2,t(x ,d)]. 



Remark 1.3. 1. Note that in the first case, t(xo,d) indeed is comparable to 
the time that F(u) has traveled from xq to xq ± de n , and thus we can say 
that the free boundary regularizes in a scale proportional to the distance it 
has traveled. 

2. Our result extends to the case where the star-shaped condition (I- a) is 
replaced by 

(I-a)' f^o is locally Lipschitz with a sufficiently small Lipschitz constant. 
We discuss the difference in the proof in this case, in section 6. 

Let us finish this section with an outline of the paper. In section 2 we 
introduce preliminary results and notations, to be used in the paper. In section 
3 we prove some properties on the evolution of solutions of (ST2) with star- 
shaped data. In addition to Harnack inequalities, we show that the solution 
stays near the star-shaped profile for a unit time (Lemma 13. 1[) . which in turn 
yields that the solution stays very close to harmonic functions (Lemma 13. 6j) . 
Making use of the results in section 3, in section 4 we perform a decomposition 
procedure to show that for a unit time all parts of the free boundary stay close 
to Lipschitz profiles, regardless of the local dynamics between the phases. This 
completes our main step in the analysis. In section 5 we describe the procedure 
leading to further rcgularization, pointing out the main difference between the 
previous results. In section 6 we discuss a generalized proof for the corresponding 
regularization result fTheoren l6.1[) when the star-shapedness of the initial data 
(I-a) is replaced by a local version (I-a)'. 

2 Preliminary lemmas and notations 

We introduce some notations. 

• For x € lR n , denote x = (x' , x n ) € -/R n_1 x IR where x n = x • e„. 

• Let B r (x) be the space ball of radius r, centered at x. 

• Let Q r := B r (0) x [— r 2 , r 2 ] be the parabolic cube and let K r := B r (0) x [— r, r] 
be the hyperbolic cube. 
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• A caloric function in n Q r will denote a nonnegative solution of the heat 
equation, vanishing along the lateral boundary of Q. 

• For xq G Tq and e n = xo/|xo|, define 

t{x ,d) := mm — — , — — — j. 

u + [Xo — ae n , OJ u (xo + ae n ,0) 

• Given e > 0, a function w is called e-monotonc in the direction t if 

u(p + At) > for any A > e. 

• W x (8 x ,e) and W t (0*, i/) with e € J?" and z/ G span(e n , e t ) respectively denote 
a spatial circular cone of aperture 28 x and axis in the direction of e, and a 
two-dimensional space-time cone in (e n ,et) plane of aperture 20* and axis in 
the direction of v. 

• w is e-monotonc in a cone of directions if w is e-monotone in every direction 
in the cone. 

• C is called an universal constant if it depends only on the dimension n and 
the regularity constant No of uq. 

The first lemma is a direct consequence of the interior Harnack inequalities 
proved in [C-C]. 

Lemma 2.1 ( [C-C] ). Suppose w(x) : lR n — > 1R has bounded Laplacian. Then 
w is Holder continuous with its constant depending on the Laplacian bound. 

Lemma 2.2 ( jFGSlj . Theorem 3). Let SI be a domain in JR n x Ft such that 
(0,0) is on its lateral boundary. Suppose Q is a Lip 1 ' 1 / 2 domain, i.e., 

n = {(x',x n ,t) : \x'\ < l,\x n \ <2L,\t\ < l,x n < f(x',t)}, 

where f satisfies \ f(x',t) — f(y',s)\ < L(\x' — y'\ + \t — s] 1 / 2 .) If u is a caloric 
function in ; then there exists C = C(n,L), where L is the Lipschitz constant 
for ; such that 

u(x,t) < u(-Le n , 1/2) 



v(x,t) ~ v(-Le n ,-l/2)' 
for (x,t) G Qi/2- 

Lemma 2.3 ([ACS1 , Theorem 1). Let Q, be a Lipschitz domain in M n x M, 
i.e., 

Qi n n = Qi n : x„ < f(x',t)}, 

where f satisfies \f(x,t) — f(y,s)\ < L(\x — y\ + \t — s\). Let u be a caloric 
function in Qi with (0, 0) G <9S1 and u(—e ni 0) = m > and supg 1 u = M . 
Then there exists a constant C , depending only on n, L, jj such that 

u(x, t + p 2 ) < Cu(x, t - p 2 ) 

for all (x, t) G Qi/2 H ^ and for < p < d x t . 
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Lemma 2.4 f |ACSlj . Lemma 5). Let u and Q. be as in Lemma [2~M Then there 
exist a, 5 > depending only on n, L, -jj such that 

w+ := u + u 1+a and W- := u — u 1+a 

are subharmonic and superharmonic, respectively, in QjflOn{i = 0}. 

Next we state several properties of harmonic functions: 

Lemma 2.5 ([D]). Let U\,u<z be two nonnegative harmonic functions in a do- 
main D of LR n of the form 

D = {(x',x n ) g -K"" 1 x 1R : \x'\ < 2,\x n \ < 2L,x n > f(x')} 

with f a Lipschitz function with constant less than L and /(0) = 0. Assume 
further that u\ = U2 — along the graph of f . Then in 

D1/2 = {\x'\ < 1, \x n \ < L,x n > f(x')} 

we have 

n , r ^ U^x'.Xn) U 2 (0,L) 

U < Oi < — r — < G 2 

U2(x',X n ) Ul{0,L) 

with G\,Ci depending only on L. 

Lemma 2.6 (iJK ). Let D, m and U2 be as in Lemma \2.5\ Assume further 
that 

m(0,£/2) = 

u 2 (0,L/2) 

Then, Ui{x' , x n )/u2(x' , x n ) is Holder continuous in D x / 2 for some coefficient a, 
both a and the C a norm ofu\/u2 depending only on L. 

Lemma 2.7 ( }C2| ). Let u be as in Lemma \2.5\ Then there exists c > depend- 
ing only on L such that for < d < c, -^—u(0, d) > and 

Cl ^<-^(0,d)<C2^ 

d ax n d 

where Q = Q(M). 

Lemma 2.8 ( |JK| . Lemma 4.1). Let f2 be Lipschitz domain contained in Biq(0). 
There exists a dimensional constant j3„ > such that for any £ € dfl, < 2r < 1 
and positive harmonic function a inflfl £> 2r (C), ^ u van i s hes continuously on 
B 2r (C) H dn, then for x G O H B r (Q, 

u(x) < C(t£)*S,up{u(y) : y € dB 2r (C) H 0} 
r 

where C depends only on the Lipchitz constants of CI. 
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Next, we point out that we use the notion of viscosity solutions for our 
investigation. When {u = 0} is of zero Lebesgue measure, it was proved 
in [KPj that the viscosity solution of (ST2) is unique and coincides with the 
usual weak solutions. (See [KPj for the definition as well as other properties of 
viscosity solutions.) Below we state important properties of viscosity solutions. 

Lemma 2.9. Suppose u is a viscosity solution of (ST2). Then 

(a) u is caloric in its positive and negative phases. 

(b) —u is also a viscosity solution of (ST2) with boundary data —g. 

(c) u + = max(u,0) (or u~ = — min(M,0) / ) is a viscosity subsolution (or su- 
persolution) of (ST2) with initial data Uq (or Uq ). 

Lemma 2.10 (Comparison principle, |KPj ). Let u,v be respectively viscosity 
sub- and supers olutions of (ST2) in D x (0, T) C Q with initial data u -< Vq in 
D. If u<v on dD and u < v on ODD Q(u) forO<t<T, then u(-,t) -< v(-,t) 
in D forte [0,T). 

Below we state a distance estimate for the free boundary and Harnack in- 
equality for the one-phase solution u of (ST1). 

Lemma 2.11 f |CKj . Lemma 2.2). Let u be given as in Theorem li.il There 
exists to = to(No, Mq, n) > such that if xq <G Tq and t <to, then 

i_£i/(2-«) < d(x ,t) < Ct 1 ^ 2 '^ (2.1) 

where a and (3 are given in (1.2]) , C depends on Nq, Mq and n, and d(xo,t) 
denotes the distance that T moved from the point xq during the time t, i.e., 

d(xo,t) := sup{ci : u(xo + de n , i) > 0}. 

Lemma 2.12 ( jCKj . Lemma 2.3). Let u be given as in Theorem li.il There 
exists do depending on No, Mo and n such that if xq € To and d < do, then 

u{xq — de n , t) < Cu(xo ~ de n , 0) for < t < t{xo, d) 

where C depends on No, Mo and n. 

The following monotonicity formula by Alt-Caffarclli-Fricdman prevents the 
scenario that both phases compete with large pressure in our problem. 

Lemma 2.13 ( |ACFj ). Let h + and h- be nonnegative continuous functions in 
-Bi(O) such that Ah± > and h+ • h— = in £>i(0). Then the functional 

r Jb t (o) \x\ J Br (o) \ x \ 

is monotone increasing in r, < r < 1. 



Corollary 2.14. Let 9Qq C M n be star-shaped with respect to i?i(0) C Slo and 
suppose #4/3(0) C Oo C i?5/3(0). Let h + be the harmonic junction in VLq — B\(0) 
with boundary values h + — on dtto, and h+ = 1 on dBi(0). Let be the 
harmonic function in £?2(0) — with boundary values h- — on dfla, and 
h— = 1 on 9£?2(0). Then there exists a sufficiently large dimensional constant 
M > such that 

b^Zl^l > M implies M*o + re„) < 1 
r r 

for S dflo, e n — x/\x\ and < r < 1/6. 
Proof. It follows from Lemma 12.131 since 

'h + (x -re n ) h^(x Q +re n ) 



1 f ' V ^' 2 dx-f , ' Vfe -,' 2 dx 



i 2r ) 4 J B r/2 (x -re n ) \X ~ X \ n 2 J B r/2 (x +re n ) \^ ~ X ^ n 

< , V^dx.f , ' Vfe -,' 2 „dx 



(2r) 4 7 B2r(:Co) |a; - x |™- 2 7 B2r(:co) |x - x Q \ n ~ 
4>{2r) < 0(1/3) < C n . 



□ 



3 Properties of solutions with star-shaped ini- 
tial data 

Lemma 3.1. Iftlo and uq are star-shaped with respect to the ball B ro (0) C fV' 
then fit(it) and u(-,t) stays a-close to star-shaped for all < t < ^cr 1 / 5 . (See 
Figure 2) 

Proof. 1. Observe that, for any a > 0, the parabolic scaling (x,t) (ax,a 2 t) 
preserves both the heat operator and the boundary motion law in (ST2). There- 
fore, for any a > the function 

ui(x, t) := u((l + cr)(x - x ) + x , (1 + cr) 2 t) 

is also a viscosity solution of (ST2) with corresponding initial data. 

2. Choose x € B ro (0). Take a small c > such that B ro+C0 (0) C f^o- We 
claim that for < S < cr 6/5 , 

ui(x,0)<u(x,S) in B R (0) - B ro+Co (0) (3.1) 
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if a is small enough. To show (|3.ip , let us introduce another function 

u(x, 0) := u((l + ^)(x- x ) + x , 0). 

Also let v~ be the solution of (ST1) with initial data Uq, and with v~ = 1 on 
dBjt(0). Then by comparison, — v~ < u and Br(0) — O t (— v~) C O t (u). Hence 
by Lemma 12.111 applied for — v~ , 

n (u) C n t (u) for < t < cr 7/6 . 

Moreover, due to our assumption, 

u(x, 0) < uq(x). 

Therefore, the maximum principle for caloric functions implies 

w(x, t) < u(x, t) 

where w solves the heat equation in the cylindrical domain D = rto(u) x [0, a 7 / 6 ] 
with initial data u(x,0) and zero boundary data on dQo(u) x [0,er 7 / 6 ]. 
Now wt solves the heat equation in D, 

w t = Aw > — C at t = 0, and w t — on dflo(u). 

Therefore we conclude that wt > —C m. D . In particular 

w(x,S)>u(x,0)-C6. (3.2) 

Next we compare u\ (x, 0) with w(x, 5). Observe that for a; € -Br(O) — B ro+CQ (0), 

ui(x,0) — u(x, 0) + j°, 2 {{x — x o) ■ Du((l + s)(x — x ) + x , 0))ds 

< u(x, 0) - C (7 7 / 6 

< u(x, 0) - Ca 6 / 5 

< w(x, S) < u(x, d) 

for < S < <7 6 / 5 , where the first inequality follows from our assumption (I- 
d) on wo, the second inequality follows if a is sufficiently small, and the third 
inequality follows from (|3.2[) . Hence we conclude (13. ip . 



3. Our goal is to prove that for < 5 < cr 6 / 5 , 

ui(x, t) < u 2 {x, t) := u{x, t + S) (3-3) 

in (B R (0) — B ra+CQ (0)) x [0, cr 1 / 5 ]. Note that the inequality holds at t = by step 
2. However, we needs a bit more arguments since we do not know yet whether 
the lateral boundary data on dB ro+Co (0) is properly ordered. 
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Suppose 

fi(ui) C fi(u) for < t < t 
and f2(ui) contacts df2(it) for the first time at i = to. Observe then that 

f(x, t) := u{x, t + S) — ui (x, t) 

solves the heat equation in J7(ui) with nonnegative boundary data for < t < to, 
with 

/(x,0)>0inS fl (0)-B ro+Oo (0). 
Indeed following the computation given above, it follows that 

/(x,0) > c a in B ro+Co (0) - B ro+ ^(0). 

On the other hand, due to the fact that wt > — C and 8 < <r 6 / 5 , we have 

f(x, 0) > (w(x, S) — w(x, 0)) + (w(x, 0) — u\(x, 
Therefore we have 

f(x,t) > on dB ro+Co (0) x [0,f ] 

if to << 1. But then this contradicts Theorem 12.101 applied to the region 
(B R (0) - B ro+Co (0)) x [0,t ]. 

4. From (|3.3p of step 3, we obtain 

u((l + a){x - x ) + x , (1 + cr) 2 t) < u(x, t + S) (3.4) 

in (Br(0) — B ro+CQ (0)) x [0, cr 1 / 5 ] for any xq € B ro (0), as long as a and S are 
sufficiently small and satisfy < 5 < ct 6 / 5 . As a result, for < t < \^ x ^ b , we 
can choose 6 = a(2 + a)t < a 6 / 5 such that 

(1 + aft = t + S. 

It follows then from (|3.4j) that the function u(-,t) is cr-monotone with respect 
to the cone of directions W x in (B R (0) - B rg+co (0)) for t <E [0, \v 1/5 \- 

( Here W x = {v G S n : v = *~ for some x G S ro (0)}.) 

F - aro| 



□ 



Remark 3.2. Observe that, due to (I — b), we have for x G To 



t(x,d) := min[ ^ f ] G [d 7 / 6 ,d 5 / 6 ] « (3.5) 

u+(a; — c(e n , (J) u (a; + ae n ,(J) 

where t(x,d) is the time it takes for the free boundary to regularize in Bd(0). 
Therefore, u{-, t) is at least d 4 -monotone with respect to W x in (Bji(0)—B ra + Co (0)) 
for < t < t[xo,d). This property will serve as a basis for our regularization 
argument in section 3. 
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Lemma 3.3. (Harnack at t = 0) Let x € To, then for all s > and for 

< t < t(x, s) we have 

u + (x — se n , t) < C\u + (x — se n , 0) 

and 

u~{x + se n , t) < C\u~{x + se„, 0) 

where e n = x/\x\ . 

Proof. Let v + solve the one-phase Stefan problem (ST1) with initial data v$ (x) = 
Uq(x). Then v + is also a solution of (ST2) with uo(x) < Vq(x), and thus by 
Theorem 12. 101 we have 

u(x, t) < v + (x, t). 

Therefore it follows from one-phase Harnack inequality applied for v + (x, t) that 

u + (x — se n ,t) < v + (x — se n ,t) < C\v + {x — se n , 0) = C\u(x — e n , 0) 

for < t < t where to = s 2 /u(x — se„, 0) > t(x, s). 

As for u~(x,t), we compare u~ with the solution v~ of (ST1) with initial 
data Vq(x) = Uq (x) and with boundary data v~ = 1 on <9£?i?(0). The rest of 
the argument is parallel to above. 

□ 

Lemma 3.4. (Backward Harnack at t = 0) Let x £ To and let s > 0. Then 
for0<t< t(x, s) 

u + (x — se n , 0) < C\u + (x — se n , t) 

and 

u~(x + se n , 0) < C\vT (x + se n , t) 

Proof. We will only show the lemma for u + . The other part follows by a parallel 
argument. Let v~ solve (ST1) with initial data and with boundary data 1 
on dBfj(Q). Then — v~~ is also a solution of (ST2) with — Vq < uo, and thus by 
Theorem 12.101 — v~ < u and 

{v~ = 0} C {u > 0}. 

Note that fl(v~) moves according to the one-phase dynamics, which has 
been studied in detail by |CK2 ;. In particular we know that £l(v~) will be 
Lipschitz at each time. Moreover, for a boundary point (x,i) £ r(i; - ) and 
d := dist(x, ro(w~)), the normal velocity V x j satisfies 

, \ i v~~ (x + 2de n ,0) ,a i £=1 .„ „. 

V x , t = \Dv-(x,t)\ !s — 2d ' < d - 1 < tFS (3.6) 

where the last inequality follows from Lemma \2. Ill Let v*(x,t) solve the heat 
equation in {v~ = 0} with initial data uo(x) and boundary data on d{v~ = 0}. 
Since 

fl{v*) = {v- = 0} C {u > 0}, 
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we have v*(x, t) <u(x,t). Moreover, for any given t > 0, v~(x, s) := v~ (yix,ts) 
satisfies the assumptions of Lemma \2. 41 Thus it follows that v~(-,t) is t a -close 
to a harmonic function in B^(x) for some a > 0, where x € IV Moreover, due 
to the assumption on the initial data, (v*)t = Ai>* > — C at t = 0. Also on 

IV), 

(v*)t/\Dv*\ = -(v-) t /\Dv-\ = -\Dv~\ > -t%£ 

where the last inequality follows from (|3.6|) . Since Sl(v*) is Lipschitz and 
r t (w*) = T t (v~) is regularized in time (Theorem ll.l[) , (|3.6p also holds for \Dv*\. 

Hence on T(v*), 

(v*) t = -|DtT||iV| > -i^ 1 > -t- 2 ' 5 . 
Since (v*)t solves a heat equation in Q(v*), it follows that for x 6 To, 

(v*) t > -t- 2/5 in B^- t/2 {x - Vte n ) x [0, t]. (3.7) 

Then since v*(x- y/te n ,0) > (Vt) a > (Vt) 7/6 = t 7/12 , 

u + (x — \/ie n , 0) = v* (x — V~te n , 0) < C\v*{x — Vie n , t) 

< C\u + {x — yie n , t) 

where the first inequality follows from (|3.7j) . Since T(v*) = r(u _ ) is Lipschitz 
in a parabolic scaling, v* is almost harmonic. Hence v*(-,t) is bigger than the 
harmonic function w*(x) in tt t (v*) n B ^q{x) with its value 

u*(x - Vte n ) = (Ci)- 1 u+(x - Vten, 0). 

Note that if < t < t(x, s), then s < y/t. Hence for < t < t(x, s), 

C\u + (x — se n , t) > C\V*{x — se n , t) > C\ui (x ~ se n ) > Cu + (x — se n , 0), 

where the last inequality follows since the one-phase result implies a power law 
on the movement of r(i>~) = T(v*) (see Lemma 2.5 of [CJK1]), and this yields 
a bound on u + (x — se„, 0)/uj t (x — se n ). 

Similar arguments apply to u~ , if we consider the function v + solving (ST1) 
with initial data Uq, and the function v* solving the heat equation in {v + = 0} 
with initial data uq and with boundary data on T(u + ) and —1 on <9Br(0). 

□ 

Lemma 3.5. (Distance estimate at t = 0) Let x € To and let s be a suffi- 
ciently small positive constant. If 

\u + (x — se n ,0)\ . , \u~(x + sen, 0)| 

< m and < m, 

s s 

then forte [0, 

d(x, t) = sup{r : x + re n or x — re n € T t (u)} < s. 
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Proof. Let v + solve (ST1) with initial data Uq , and let v~ solve (ST1) with 
initial data Uq and with v~ = 1 on dBn(0). Then by comparison, — v~ <u< 
v + and the lemma follows from the one- phase result Theorem ll.il □ 

In the following lemma, we approximate our solution by harmonic functions. 

Lemma 3.6. (Spatial regularity in the whole domain) For xq £ r and 

r > 0, there exists a function uj(x,t) := u) + (x, t) — U)~(x,t) such that 

(a) ui(-,t) is harmonic in its positive and negative phase in 

(1 + r)Q t (u) — (1 — r)Q t (u), and Q,(lu + ), Q(oj~) are star-shaped; 

(b) For a dimensional constant C > 0, we have 

u + (x,t) < u+(x,t) < Coj+{{1 - r 5 ^)x,t) 

and 

(j-(x,t) < u~(x,t) < Cuj-((l+r 5/4 )x,t) 

in B r (x ) x [r 2 ,t(x ,r)}. 

Remark 3.7. Note that we do not know yet whether the solution is close to a 
Lipschitz graph in time. Also, note that t(xo,r) > ?* 7 / 6 3> r 2 , and d{u} + > 0} 
need not be 3{uj~ > 0}. 

Proof. 1. We will only show the lemma for u + . Let T* be the free boundary 
obtained from the one-phase problem (ST1) with the initial data Uq", and let 
O* be the region bounded by T*. Let V\ solve the heat equation in fT and in 
-Br(O) x [0, 1] — Q*, with initial data uq and with v\ = —1 on dBn(0). Similarly, 
we define Vi , whose free boundary is obtained from the one- phase solution with 
initial data . Then by comparison, 

V2 < U < Vl- 

Hence the free boundary of u is trapped between the free boundaries of v\ 
and V2- Also, since one-phase versions v\ and v-i behave nicely, we have those 
functions almost harmonic up to r-ncighborhood of their free boundaries for 
r 2 /2 < t < r 2 . Next note that the range of t is < t < t(xo,r), and thus 
both of the sets T t (vi) and T t (v2) are within distance r of To(u) in B t {xq). In 
particular, using the one-phase result, i.e., arguing as in Lemmas 2.1 and 2.3 of 
[CK], we obtain 

v 2 {x - 2re n ,t) - u (x - 2re„,0) - v\(x - 2re n ,t) (3.8) 
for < t < t(x ,r). 
2. Observe that 

t(x ,r) < r 2 ~ a < r 5/6 := r. 

Due to Lemma [3~Tl we know that at each time, U t (u) is t 5 -close to a star-shaped 
domain D t up to the time t = r, i.e., 

D t C n t (u) C (1 + r 5 )A C (1 + r 4 )A (3.9) 
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Figure 2: Approximation of the positive phase by a star-shaped domain 

for < t < t. 

Then by Lemma EH1 and (|3"l?]l with /3 > 5/6, 

u(ar, t) < r^ 20 ^ = r 13 / 24 on 9(1 - r 13 ^°)D 

for < t < t. (Here note that we can apply Lemma 13.31 up to the time r since 

t(£,r 13/20 ) > r 13 ( 2 -«/ 20 > r for any z G r .) 

Then by the T 5 -monotonicity of u, 

u(x, t) < r 13 / 24 on B R (0) - (1 - r 13 / 20 + r 4 )Z? (3.10) 

for < t < t. Since T t (u) stays in the r 5 / 6 -neighborhood of Tq(u) up to r, 
we obtain that 9Dj stays in the r 25 / 36 -neighborhood of dDo up to the time r. 
Since r 25 / 36 < r 13 / 20 , (pOO]) implies 

u(x,t) < r 13/24 on B R (0) - D s (3.11) 

for any < s, t < t. 
3. Let 

to = < <i = r 2 < t 2 = 2r 2 < ... < t ko = k r 2 < t 
and fix a number b such that 

5/4 < b < 61/48. 
We will construct a supersolution of (ST2) in 

(B R (0)-{l+r b )D tk ) x [t k ,t k+1 ]. 
Let w k (x) be the harmonic function in 

E:= (l + 4r b )A fc -A fc 
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with boundary data zero on <9(1 + Ar b )D tk and C„r 13 / 24 on dD tk , where C n is 
a sufficiently large dimensional constant. Extend w(x) = in JR n — £. Next 
define 

r b-2 

${x, t) := uu>(y) :\x-y\< r b - (t - **)— } 

in (Bfl(O) - (l + r b )A fc ) x [i fc ,i fe+1 ]. We claim that u is a supersolution of (ST2) 
since our constant b satisfies 

r 6 - 2 > r !t-i>. (3.12) 

To check this, first note that $(-,£) is superharmonic in its positive set and 
$ t > 0. Hence we only need to show that 

^ > |£$| on (3.13) 

Due to the definition of T t ($) has an interior ball of radius at least r b /2 for 
tk < t < ife+i. This and the supcrharmonicity of $ in the positive set yields 
that 

£ 13/24 

\D$\ < — - b — on T($) 

for a dimensional constant C > 0. Moreover T($) evolves with normal velocity 
\r h - 2 . Since ([3~T2]) holds for our choice of b (i.e., for 5/4 < b < 61/48), we 
conclude (|3 . 1 3p for r smaller than a dimensional constant r(n). Now we compare 
u with $ in 

(B R (0)-(l + r b )D tk ) x [t fc ,f fc+ i]. 

Note that by (|3~TT|) , 

u + < $ on 9(1 + r b )D tk 
if C„ is chosen sufficiently large. Also at t = tk, (|3.9|) implies 

«(-,tfc) < < on B R (0) - (l+r b )D tk . 

Hence we get u < $ in (JR™ - (1 + r b )D tk ) x [t t ,t t+ i]. This implies 

C n($) U ((1 + r b )D tk x [t fcj tfc+i]) := fl($) (3.14) 

for t k <t < tk+i- 

4. Next we let v(a;, solve the heat equation in 

n($)-((l-3r)n (u) x [t fc ,t fc+ i]) 

with initial data v(-,tk) — u(-,tk) and boundary data zero on T($) and v — u 
on (1 — 3r)To(u). Observe that, due to (|3.14p . we have u + < v for tk < t < tk+i- 

Since is star-shaped and expands with its normal velocity < r b ~ 2 which 
is less than r^ 1 , Lemma 12.41 applies to v(x,t) := v{rx,r 2 t). In particular there 
exists a constant C > such that 

{1/C)v(x,t) < hi{x,t) < Cv(x,t) 
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for (t k + ifc+i)/2 < t < tk+i, where hi(-,t) is the harmonic function in 
flt(v) — (1 — 2r)f2o(w) with boundary data zero on T t (v) and v on (1 — 2r)To(w). 
Hence we conclude that 

u + < v < Chi 
in (S fl (0)-(l-2r)no(«)) x [(t k + t k+1 )/2,t k+1 }. 

5. Similar arguments, now pushing the boundary purely by the minus phase 
given by the harmonic function yields that 

B H (o)-n t (*) :=n t cfit(«) 

for tfc < t < tfe+i, where 

n t = {x£ D tk : dist(x, 9AJ > 3r b + — (i - i fc )}. 
Let w(x,t) solve the heat equation in 

H - ((1 - 3r)Jlo(tt) x [tk,t k +i})) 

with initial data ife) and boundary data zero on <9IT, and uon (1- 3r)To(u). 
Then u > t). 

Since IT is star-shaped and shrinks with its normal velocity < r b ~ 2 which is 
less than r" 1 , Lemma T2.4I applies to w(x,t) := w(rx,r 2 t). In particular there 
exists C > such that 

u+ > w > (1/C)h 2 

for (tj. + tk+i)/2 <t< ife+i, where h,2(-,t) is the harmonic function in 
II f — (1 — 2r)flo(u) with boundary data coinciding with that of w. 

6. Lastly we will show that hi and h 2 are not too far away, i.e. 

hi(x,t) < Ch 2 (x-8r b e„,t) 

with a dimensional constant C > 0. Since u is between (1/C)h 2 and CTii, 
this will conclude our lemma for (tk + t k +i)/2 < t < t k +i- Then by changing 
the time intervals [tk,t k +i] to [tk + r 2 /2,t k +i + r 2 /2], we obtain lemma for 
r 2 < r < t(xo, r). 

To prove the claim, observe that 

n t ( w ) c n t (v) c (i + 8r 6 )OtH 

Moreover, observe that 

v 2 ((l + r b )x,(l + r b ) 2 (t-t k )+t k ) < v(x,t),w(x,t) 

< Vl ((l-r b )x,(l-r b ) 2 (t-t k )+t k ) 

for t k <t< t k+1 . This and (O yield 

v(x - 2re n ,t) ~ w(a;o - 2re n ,t) ~ m(x - 2re„,0). 



18 



It follows that 

w(x,t) < v(x,t) < Cw(x - 8r b e n ,t) on (1 - 2r)T x [t k ,t k+1 ]. 

Hence due to Dahlbcrg's lemma, we conclude that 

hi(x,t) < dv{x,t) < C 2 w{x - 8r b e ni t) < C 3 h 2 (x - 8r b e n ,t) 

in B r (xo) x [(t k + tfc+i)/2, tfc+i]. Since the inequality holds for any 5/4 < b < 
61/48, we can conclude the lemma. 

□ 

Proposition 3.8. (Regularization in bad balls) For a fixed xq G To(u), 
suppose that either 

u + {xq - re n ,t ) > Mu~(x + re n ,t ) 

or 

u~(x + re„,t ) > Mu + (x a - re n ,t ) 

for M > M„, where M n is a sufficiently large dimensional constant. Then for 
r < l/M n) there exists a dimensional constant C > such that 

IT? +i 4-w s- r< u+ ( x ° ~ re «'*o) , [V7 _/ ,x, . n u~{x +re n ,t ) 

Ywv/iXjt) \<C and \vu (x,t)\<C 

r r 

in B r (x ) x [t(x ,r)/2,t(x ,r)]. 

Remark 3.9. 1. In the next section, we will extend this Lemma for later times, 
i.e., for xo G T tg . (See Lemma \%71\ ) 

2. Note that the situation given in Proposition \3.8\ is essentially a perturba- 
tion of the one-phase case in \CKf . The main step in the proof is in verification 
of this observation: i.e., by barrier arguments we will show that our solution is 
very close to a re-scaled version of the one-phase solution of (ST), for which the 
regularity of solutions are well-understood (see Theorem 1.1). 

Proof. Without loss of generality, we may assume that 

u + (xq — re n , 0) > Mu~ (x + re n , 0) . 

1. First we show that after a small amount of time u become almost harmonic 
near the free bounadry. By Lemmas I3.3l and I3.4l implv that for < t < t(xQ,r), 

u + (xq - re„,t) — u + (xq - re n , 0), u"(x Q + re n ,t) ~ u~(x Q + re„, 0) (3.15) 

Also note that, by the assumption on the initial data Uq, Lemma 13.61 holds at 
t = 0. In other words, there exists a function u)(x, 0) = u)q{x) such that 

(a) luq is harmonic in its positive and negative phases in 
(l + r)no(u)-(l-r)no(u); 
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(b) r2(wg~) and Q(ui ) are star-shaped; 

(c) In B r (xo) , we have 

wjj"(a?) < < Oj+((1 - r 5 / 4 )a;) (3.16) 

and 

Ua(x) < u-(x) < Cu:^((l + r 5 ^)x). (3.17) 

Next we improve p,16|) and (|3.17p for later times, and obtained the inequal- 
ities with C = (1 + r a ) for t > r 3 / 2 . By the distance estimate-Lemma 12 . 1 1[ the 
free boundary of u moves less that r 9 / 7 < r 5 / 4 during the time t = r 3 / 2 . Then 
we let v i solve the heat equation in cylindrical domains 

(l + 2r 5 /4)f] ( w +) x [0,r 3 / 2 ]U(B 2 (0)-(l + 2r 5 / 4 )tt (^ + )) X [0,r 3/2 ], 

with initial data uq and lateral boundary data zero on (1 + 2r 5 / 4 )r (w+) x 
[0, r 3 / 2 ], and —1 on dB 2 (0) x [0, r 3 / 2 ]. Similarly, we let v 2 solve the heat equation 
in cylindrical domains 

(1 - 2r 5 / 4 )ft (u; + ) x [0,r 3 / 2 ] U (B 2 (0) - (1 - 2r 5 / 4 )ft (w+)) x [0,r 3 / 2 ] 

with initial data uq and lateral boundary data zero on (1 — 2r 5 ' 4 )ro(w + ) x 
[0, r 3 / 2 ], and —1 on dB 2 (0) x [0,r 3 / 2 ]. Then by comparison, w 2 < u < V\. Also 
by Lemma O and ^ > 5/6, 

|«i-«a| <^ fxf =r 25 / 24 
in the domain. Note that on (1 — r 6 / 7 )ro(w + ), |ui| > r? x ^ = r and thus 
\vi -v 2 \< r ai \ Vl \ on (1 - r 6/7 )r (u; + ) for a x = 1/24. 

Similarly, 

\vi-v 2 \ < r ai \v 2 \ on (1 + r 6 / 7 )r (w+). 

Then since v\ and v 2 are almost harmonic in the r 3 / 4 -neighborhood of their 
boundaries for ^r 3 / 2 < t < r 3 / 2 , the above inequalities on \vi — v 2 \ imply the 
following: for ^r 3 / 2 <t< r 3 / 2 , there exist positive harmonic functions £j + (-,t) 
and u)~(- 7 t) defined respectively in 

n t (v+) n (B R {0) - (1 - r 1 - 6 )O (^ + )) and n t («7/) n (1 + r 1 " 6 )^ (w + )) 

where b = 1/7, such that for some a > 

t) < u+(x, t) < (1 + r a )w+((l - 4r 5 / 4 )x, t) (3.18) 

and 

uT(x,i) < u~(x,t) < (l + r a )u-((l+4r 5/4 )x,t). (3.19) 
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Now on the time interval [0,r 3 / 2 ] + |r 3 / 2 , 1 < k < m, we construct v\ and 
x>2 so that they solve the heat equation in the cylindrical domains with 

r(i*) = (1 + 2r 5 / 4 )r |l , 3/2 (- + ) x [|r 3 / 2 , (1 + k -)r^] 

and 

r(« 2 ) = (1 - 2^)1^3/^+) x [^r 3 / 2 , (1 + ^)r 3 / 2 ]. 

Then by a similar argument as above, we obtain harmonic functions Co {•,£) 
satisfying (|3~18| and (f37T9|) for 

ii-^r 3 / 2 <t< (l + ^)r 3 / 2 . 

Hence we conclude (j3~Tg|) and ([3~Ti?j) for r 3 / 2 < t < t(x 0l r). 
2. Next we re-scale u(x,t) as follows: 

u(x,t) :— a~ 1 u(rx + Xo,r a~ t) in 2Q Xo , 
where a := u + (xq — re„,<o) < 1- Then u(x,t) solves 

' (ad t - A)u = in Q(u) 

V = \Du + \ - \Du~\ on T{u) 

< 

u(-e„,0) = 1 

u(e n ,0) = -l/N where N > M. 
Furthermore, (|3.15|) implies that for < t < 1, 

£t + (-e„,i) ~ 1, u~(e n ,t) ~ — . 

Let u> be the corresponding re-scaled version of u> given in f|3 . 1 8[) and (|3.19[) . 
then in B r -b(0) n f2o(w) we have 

(1 - r a )w + ((l +4,r 5/4 )x,ar- 1/2 ) < i + (i,ar" 1/2 ) < w + {x,ar~ 1/2 ) (3.20) 

and 

(1 -r")r(.i,ar 1/2 ) <r(a;,ar 1/2 ) < uT((i + 4r 5/4 );r, «r 1/2 ) (3.21) 
Here note that 

-1/2 /- u+ ( x o ~ re n ,t ) 1/3 
r 

Lastly, for given xq <E T(u) D -Bi(O), a similar argument as in (|3.7[) implies that 

u(x,t) < (l + r b )u(x,0) in dBi r - b (r- b e n ) x [0,1]. (3.22) 
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3. We claim that we can construct a supersolution U\ and a subsolution U2 
of (ST2) such that 

U 2 (x,t) < u(x,t) < Ui(x,t) < U 2 (x- y/een,t) in Bi(0) x [ar- 1/2 ,l] 

and that U2 is a smooth solution with uniformly Lipschitz boundary in space 
and time. Then for sufficiently small r > the lemma will follow from analysis 
parallel to that of |ACS2j . 

To illustrate the main ideas, let us first assume that 

(a) (|3.20p and (|3.21[) hold in the entire ring domain R x [0, 1], where 

R={x:d(x,T Q (u)) <r- b }; 

(b) u(x, t) <{1 + r b )u(x, 0) on dR x [0, 1]. 

Let E/+ be the solution of (HS) in £ = (JR n - (O - R)) X [0, 1] with initial 
data w(x, t) and boundary data (1 + r b )u(x, 0), and let 

Ut = - C7f in R x [0,1], 

where (■, t) is the harmonic function in R— J7( with fixed boundary data 
zero on r([/+) and C/N on dR - r2(C/ 1 + ). Then U x is a supersolution of (ST2) 
in S, and thus by Theorem 12.101 and the assumptions (a)-(b) we have u < U\ 
in E. 

4. The construction of the subsolution U2 is a bit less straightforward. We 

use 

U+(x,t) :=(l-e) sup C/+((l + y^Jy, t), 

|K-x|<V?(l-c(i)) 

where e=l/N and c(*) := i 4/5 . Then we define 

where i? is the ring domain as given above and U2 (•, t) is the harmonic function 
in R— £l(U2 ) w hh fixed boundary data zero on r(/7^~) and C/N on dR— 
Then [/2 satisfies the free boundary condition 

Vfr a < (l + e)|^C/ 2 + |-V^c'(i)- 
Therefore, J7 2 is a subsolution of (ST2) if we can show that 

Vec'(t) > e\DU^\ + \DU 2 ~\ on T{U 2 ) (3.23) 

and Jg 1 c'(s)ds < 1. 

The analysis performed in CK , as in the proof of (c) of Theorem ll.il yields 
the following: at a fixed time t, T(Ui) regularizes in the scale of d := d(t) which 
solves 

d 2 

*~ f/i(-de n ,0)' 
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T lie re fore 

Ml ~ Elt^l and \ DU -\ „ 

on 

r(c/ 2 ) x [t/2,t]. 

Observe that since /3 > 5/6, 

Uf(-de n ,0) < d 5/e and JXf (de n ,0) < ed 5/6 , 

then we have 

^ t/ 2 + (-de n ,0) + £/ 2 ~(de»,0) < ^_ 1/6 < /^-i/ 5 
rf d 
where the last inequality follows from 

t = d 2 /Ui{-de n ,0) < d 2 /d a < d 5 ^. 

Hence c{t) — i 4 / 5 satisfies (|3.23p . and we conclude that U 2 is a subsolution of 
(ST2). 

Now we can use the fact 

U 2 < u < Ui in B c (0) x [0,c] 

to conclude that u is y 7 ?- close to U\: a Lipschitz (and smooth) solution in 
Bi(0) x [1/2, 1]. Once we can confirm this, everything else follows from analysis 
parallel to that of [ACS2] with the choice of a sufficiently small e. 

5. Now we proceed to the general proof without the simplified assumptions 
(a) and (b) in step 3, which are replaced with local inequalities (|3.20[) - (|3.21[) 
and (|3.22l) . For this we need to perturb the initial data outside of -Bi(O) (see 
section 4, p 2781-2783 of [CJK2] ). to obtain functions W x (x) and W 2 (x) which 
satisfies the followings: 

(a) {Wk > 0} with fc = 1,2 is star-shaped and coincides with Q ar -i/i(w) in 
B r -»(0); 

(b) {W 2 > 0} C n ar - 1/2 (w) C {W! > 0} ; 

(c) d(x, {W k > 0}) > r- b with k = 1, 2 for x € T ar -y 2 (w) n (M n - B 2r - b (0)); 

(d) Wk is harmonic in {Wk > 0} — K with boundary data zero on T{Wk) and 
(1 + r b )w{x, ar -1 / 2 ) on where 

K = {x:d(x,T(Wk))>r- b }. 

Let C/fc be the solution of Hele-Shaw problem in 

M n -^{W k >0}x[ar- 1 / 2 ,1] 

with initial data W\ and with lateral boundary data (1 + r h )w(x, ar^ 1 / 2 ). Due 
to Proposition 4.1 of jCJK2j . for sufficiently small r > 0, the level sets of Ui is 
then ec-close to those of U 2 in -Bi(O) x [0, 1]. Hence we can use U 2 instead of Ui 
in step 4. and proceed as in step 4 to conclude. □ 



23 



4 Decomposition based on local phase dynamics 

Throughout the rest of the paper, we fix xq £ Tq and a sufficiently small C( 
stant r > 0, and will prove the regularization of the solution in B r (xo) 

\+(rr. r-W Wo alar, flv a rrmafant M > M wnoro M ic a an 



con- 

x 



sranr r > u, ana win prove ine reguiarization 01 me soiuuon in n T \x§) >< 
[t(xo,r)/2,t(xo,r)]. We also fix a constant M > M n , where M n is a suffi 
ciently large dimensional constant. If the ratio between u + (xq — re„,0) and 
u~(xq + re n , 0) is bigger than M, then we can directly apply Proposition 13. 81 to 
prove the main theorem. Therefore we assume that 

Ar 1 u~(x + re„,0) < u + (x - re„,0) < Mu~(x + re„,0). (4.1) 

Let 

u + (x - re n ,0) u~(x + re n ,0) 

Co := max[ , J. 

r r 

Then since Uq and Uq are comparable with harmonic functions, Co is less than 
a constant depending on n and M (See Corollary 12 .14[) . Also note that 

C >r a_1 > r 1/6 . 

Let 

A+ = {x e T n S 2r (a; ) : u+ ( x ~ se "'°) > M( j Q for gomc r 5/4 < s < r } 

s 

and 

A-={ier„n £ 2r (z ) : - ^ + se "'°) > M C for some r 5 / 4 < s < r}. 

s 

Denote 

Lemma 4.1. 7f 

M ± (a; =p se„, 0) 

> A/Co /or some s < r, 

s 

i/iera 

u =F (a;± se„,0) 
s 

Proof. Since are comparable with harmonic functions h^- , we can argue 
similarly as in Corollary 12 .141 Observe 



;£(x - se n ) u (x + se n ) ^ h + (x - se n ) h (x + se n ) 
s s s s 



□ 
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Figure 3: Decomposition of the domain 



Now for x E A + , we can find the largest constant r x < r such that 

" +( "* e "' 0) = A/Co 
r x 

then let 

Also for x G A~ . we can similarly define r x and Q x . Let 



E :=-B r (x ) x [0,t(.T O ,r)] - |J Q a . (4.2) 



(See Figure 3) 

The following statement is a direct consequence of the definition (|4.2D . 
Lemma 4.2. J/x e T n E 0; tten /or a/1 r 5/4 < s < r 
u+(x - se n ,0) u~(x + se„,0) 

j S MUq. 

s s 

The next proposition is the main result in this section, which states that the 
solution is "well-behaved" in E. 

Proposition 4.3. There exists a dimensional constant K > such that for all 

{x,t) e rn E 

f a\ u+(x~se ni t) u-(x + se n ,t) 5/4 

(A) , < KM Go for r ' < s < r. 



Before proving Proposition 14.31 we show an immediate consequence of the 
proposition: we are ready to show that T(u) is close to a Lipschitz graph in time 
as well as in space. 
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Corollary 4.4. for (x, tjeTflS, suppose (x + ke n , t + r) G T. Then 

5/4 

' fc '^ 5/4 ^^i^" 

where K\ is a dimensional constant. 

Proof. Due to Lemma l3.6[ at any time < t < i(xo,r), we have 

/^(rc.t) < i^fof) < Ci/i^iT^e^f). (4.3) 

in B r (x$), where h := h + (-,t) — h~(-,t) is harmonic in its positive and negative 
phase in (1 + r)Clt(u) — (1 — r)fl t (u), and the domains tt(h + ) and 0(/i~) arc 
both star-shaped with respect to B ro (0). Let us pick (yo,to) G L n £. Due to 
Proposition 14. 3[ (|4.3p and the Harnack inequality for harmonic functions, we 
have 

sup u(y,*o) < CdKMC r 5 / 4 (4.4) 
where C is a dimensional constant. On the other hand, due to Lemma 13. II and 

£5 < r 25/6^ we naye 



«(-,*o) < in Bi r 5/4(y + r 5 / 4 e n ) 



(4.5) 



Let 



2/1 := 2/0 



+ r 5 / 4 e„, C 2 := CC X KMC^ r(t) := ^r 5 / 4 - C 3 (f - to) 



where C3 = CC2. Next we define 4>(x,t) in the domain 

r 5/4 

II := B 2r5/i ( yi ) x [t , t + -7H 

03 



such that 



-A0(-,t) = in B 2 r 5 / 4 (yi) - B r(t)(Vi) 
= 2C 2 r 5 /4 n aS 2r5/4 ( yi ) 



,0 = in B r{t) (yi). 

Then by (|P|) , and (|4"3|) . it -< </> at t = t in LT. Let T Q be the first time 

where u hits </> from below in LT. Since (|4.4p also holds for any (x,t) G L n £ in 
place of (yo, t ), we have u < <j) on the parabolic boundary of Lin {to < t < To}. 
On the other hand, if C is chosen sufficiently large, then 

6t r 5/4 
7^7 - C 3 > \D<f>\ on dB r{t) { Vl ) x [t 0j *i := to + ^r], 

and thus ^ is a supersolution of (ST) . This and Theorem 12.101 applied to u 
and <fi in LT yields a contradiction, and we conclude that L(w) lies outside of 
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Bi r 5/4.(ya + r 5 ' 4 e„) for to < t < t\. Similarly, by constructing a negative 
radial barrier and comparing it with u, one can show that T(u) lies outside of 
Bi r s/i(yo — r 5 / 4 e n ) for io <t<t%. Hence we conclude. 

□ 

We proceed to show our main result, Proposition ^. 31 The following lemmas 
are used in the proof of the proposition. 

• For xq <E T to , define 

t{x$,r) := min| 



u+(x - re„,t )' u {x +re n ,t )' 

Lemma 4.5 (Harnack at later times). Fix s <G [r 5 / 4 ,r]. // (yo,to) <G T n S, 
then 

u + (y - se n ,t ) > ciu + (y Q - se„,t + r) 

and 

u~(?/o + se„,t ) > ciu _ (y + se„,i + r) 
/or < t < t(j/ , s)/2 and ci > 0. 

Proof. We will show the lemma for u + : the statement on u~ follows via parallel 
arguments. 

1. Let {y ,t ) ernE and let s £ [r 5 / 4 ,r]. Let h+ be given as in (14.31) . Due 
to Lemma 13.31 and Lemma |3.4[ we have 



h + (y - 2re n ,h) < u + (y - 2re„, ti) 

< Cu + (y -2re n ,t 2 )<Ch + (y Q ~(2r + r 5 / 4 )e n ,t 2 ) 

for < t\,t<i < t + t(y ,r)/2. (Here note that yo € B r (xo).) In particular 

U +(y - 2re n ,t) < Ch + {y Q - (2r + r 5 / 4 )e„, i ) < dh + (y Q - 2re n ,t ) (4.6) 

fort <t + t(y 0l s)/2. 

2. Now let v + solve (ST1) in (M n - (1 - 2r)A ) X [*o,*o + t(y ,s)/2] with 
initial and boundary data C2h + (x — 2se n ,t). Since s > r 5 / 4 , (|4.3I) implies 

flt(u) C n t0 (v + ) C n t (v + ) in B 2s (2/o) X [t ,*o + %o,s)/2]. (4.7) 

Then by l|47f|) . (|L6|) and (|L3]) . 

u+ < i;+ in 5 s (y ) x [t , to + t(y , s)/2] 
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if we choose Ci as a multiple of C\ by a dimensional constant. Moreover, due 
to the Harnack inequality for one-phase (ST1), one can conclude that 

u + (yo - se„,i + r) < v + (y - se„,t + r) 

< Cu + (?/o - se„,t ) 

= CC 2 h + (y Q - 3se„,i ) 

< C 3 h + (y - se n ,t ) 

< C 3 u + (y - se n ,t ) 

s 2 

for < t < — — ~ t(yo,s)/2. Here the first inequality uses 

v+(ya - se n ,to) 

u + < v + , the second uses the Harnack inequality for v + , the third one uses 
the Harnack inequality for harmonic functions and the last one uses (|4.3[) . □ 



Lemma 4.6 (Backward harnack). Suppose that (A) holds up to time t = 
To < t(x ,r). If (yo, t ) G r and t < T , then for < r < t{y , s)/2, 

u + (y - se n ,t ) < Cu + (y Q - se n ,t + r) 

and 

u~(yo + se n ,t ) < Cu~(y + se n ,t + r) 
where < s < r and C is a universal constant. 

Proof. We will show the argument for u + , due to the symmetric nature of the 
claim. The argument here will be similar to that of Lemma 13. 4| replacing the 
initial data Uq and Uq (used in the construction of barriers) by h + (x,to) and 
h~(x,t ) given in (|4.3[) . 

We consider v\. a one-phase solution of (ST1) in 

n := (l + r)Q to x [t ,t +t(y ,s)/2] 

with initial and lateral boundary data C\h~ . Then v± < u in n. Now let vi 
solve the heat equation in \v\ — 0} x [to, to + t(yg, s)/2] with initial data 



«2(-,*o) = 




{ U i(-,io) = 0}-(l-r){/i+(-,i )>0} 
(l-r){fc+(-,tb)>0}, 



where ft-(-) is a C 2 extension function of h + (-, to) chosen so that h(-) < u + (-,to). 
The rest of the proof is the same as that of Lemma 13.41 

□ 
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Lemma 4.7. (Regularization in bad balls) For a fixed (xo,to) G and 
suppose 

u + (x - re n ,t ) > Mu~(x + re n ,t ) 

or 

u~(x + re n ,t ) > Mu + (x Q - re n ,t ) 

for M > M n , where M n is a dimensional constant. Then for r < 1/M„, there 
exists a dimensional constant C > such that 

ivy +, ^ n u + {x Q - re n ,t ) _ u-(x +re n ,t ) 

Vu < C and \vu \ < C 

r r 

in B r (x ) x [t + t(xo, r)/2, to + t(x , r)]. 

Proof. The proof of this lemma is parallel to that of Proposition 13.81 We use 
Harnack and backward Harnack inequalities (Lemmas 14.51 and 14. 6p instead of 
Lemmas 13.31 and 13.41 Also we have Lemma 13.61 □ 



We are now ready to prove our main result, Proposition 14.31 Observe that 
(A) holds up to some To > by Lemma 14.21 and Lemma 13.31 



Proof of Proposition \4-S\ Let K be a sufficiently large dimensional constant 
such that K ^> M. Let us assume that (A) breaks down for u + for the first 
time at t — To- Then 

" +( "°- Se "' ro) = KM C (4.8) 
s 

for some (jzo, T ) G T n S and r 5 / 4 <s<r. Let 

h = sup f h : u+ (zo-ke n ,T ) ^ ^ ^ g ^ fc ^ 
k 

Note that h < r/2 due to Lemma T3. 31 and the definition of Co, and h > 2s due 
to Lemma 13.61 By the definition of h we have 

h 

Let us find to- the closest time before To such that for some (yo,to) G T 

To - 1 = t(y , h)/2 and 2/o/|s/o| = z /\z \. 
Then Lemma |4~51 implies 



u + (yo - he n ,t ) ^ u + (y - he n ,T ) ^ u + (z - he n ,T ) _ 2 

h h h 

Since u + (-, to) and u~ to) are comparable to harmonic functions (Lemma [33 
a similar argument as in Lemma 14 . 1 1 implies that 



u (y + he n ,t ) < < 1 u + (yp - he n ,t ) 
h ~ ~ M 2 h 
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Hence by Lemma T4.71 we have 

|V U + (-,T )|~A/ 2 Coin B h (y ) 
Since B s (z ) C Bh(yo), this would contradict (|4.8p since K S> M. 



□ 



5 Regularization after t = t(r). 

Recall that xo £ To and r > are fixed, and they satisfy (|4.ip . Our goal is to 
prove the regularization of the free boundary after the time t(xo, r)/2 in B r (xo). 
Define 

Q r (x ) := B r (x Q ) x [t(x Q ,r)/2,t(x Q ,r)] C E. 

Let us briefly review the information we have on u so far. Due to Lemma 3.6 
and Corollary 4.4, our solution u is e- monotone in Q r {xo), with respect to a 
space and time cone, where the space cone W x (e n , Oq) satisfies 

|0o-7t| =0(L), 

where L is the Lipschitz constant of the initial domain Hq given by (jl.ll) . More- 
over, due to Proposition 4.3, u does not grow too big over time, which along 
with Lemma 3.8 guarantees that there is no big flux of u coming in from outside 
of B r (xo) to perturb our solution. Therefore the theory developed in jACSlj - 
|ACS2j . which says localized solutions with flat free bondaries are smooth, ap- 
plies with appropriate modifications if we have L small enough such that the 
waiting time phenomena as seen in [CK2] is prevented. More precise description 
of the situation as well as precise modifications are detailed below. 

As a result of Proposition 14.31 (A) holds up to 

t = t(x ,r) <Cr 2 - a <r 3 / 4 . 

Moreover Q r (xo) C E, and thus Corollary |4.4l and Lemma l3~Tl the free boundary 
r(it) is r 4 / 3 -monotone in Q r (xo) with respect to the time cone 
W t (e n , tan _1 (l/iir 1 MCo)) and the space cone W x (e n ,9 ). Here #o is the angle 
corresponding to the Lipschitz constant of Po, and t(xo,r) = -^r. 

On the other hand, by Lemma 13.31 and the definition of Co, 

/ t(xQ,r)\ 

u{xo - re n , y 2 ) 

Since Q r (x ) C E, Proposition 14.31 implies 

^0- < KM in B r (x ) x [t(x ,r)/2,t(x ,r)}. 



30 



Motivated from the above estimates, we consider the re-scaled function 

u(x,t) := ——u{rx + x ,r t H ). 

C r 2 

The main difficulty in applying the Method of [ACS]-[ACS2] lies in the fact 
that we cannot guarantee the e-monotonicity of the solution u in time variable 
(although we can obtain, as above, the r 4 / 3 -monotonicity of the free boundary 
T(ii)). In [ACS]-[ACS2], it was important that initially the time derivative of 
the solution was assumed to be controlled by the spatial derivative, i.e., 

\u t \<C(\Du + \ + \D U -\). (5.1) 

Using (|5.1[) one can prove that the direction vectors 

Du + ( Du- n 

do not change much for < t < I. This is pivotal in regularization procedure 
since then T(u) regularizes along the direction of the "common gain" obtained by 
those two direction vectors, the regularity of T(u) then makes above two vectors 
line up better in a smaller scale, which contributes to further regularization of 
T(u) in a finer scale. In our case we do not have (|5.1|) . which requires an extra 
care in showing that the vectors do not change their directions too rapidly. 

o e-monotonicity ofT(u) to full monotonicity of u 

First we prove that the e-monotonicity of T(u) improves to Lipschitz conti- 
nuity. Let a = Cor. Then in the domain -Bi(O) x [— ~, -], u(x,t) solves 

u t - Au = in {u > 0} 

V = a(\Du+\-\Du-\) on d{u > 0}. 

Here note that r 7 / 6 < r a < a < r^ < r 5 / 6 . In this scale, since u is Caloric and 
r(u) is r 1//3 -close to a Lipschitz graph in space and time, it follows that so does 
uin B 1/2 (0) x + 

Note that in above step we are losing a lot of information over time: T(u) is 
in fact r^-close to a Lipschitz graph moving very slow in time, but this does 
not guarantee that u also changes slowly in time. 

We then follow the iteration process in Lemma 7.2 of [ACS] to show the 
following: 

Lemma 5.1. If r is sufficiently small, then there exists < c, d < 1/2 such that 
the following is true: u is Ar 1 / 3 -monotone in the cone of directions W X {9 X — 
r d ,e„) andW t (6 t -r d ,v) in the domain Bi_ r o(0) x [ ^^- , \\. 

One can then iterate above lemma to improve the e-monotonicity to full 
monotonicity, and state the result in terms of u: 
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Lemma 5.2. u is fully monotone in #1/2(0) x [0, ^] for the cone 
Ci := W X (9 X - r d , en) U W t {9 t - r d , v), 



for some constant < d < 1/2. 
o Further regularity in space 

Now we suppose u is Lipschitz in space and time. Then in particular, we 
have the Lipschitz regularity of u in space (and very weak Lipschitz regularity 
of u in time.) We are interested in proving the following type of statement: 

Lemma 5.3 (enlargement for the cone of monotonicity). There exists A > 
such that the following holds: Suppose u is Lipschitz with respect to the cone of 
monotonicity A x (e n ,6o) in -Bi(O) x [— -, —]. Then in the half domain -81/2(0) x 
[ — 5a' Y~1> " * s lipschitz with respect to the cone of monotonicity A x (v, (1 + X)6q) 
with some unit vector v. 

To prove the enlargement of the cone, we take a closer look at the change of 
u over time, in the interior region. More precisely, we need the following lemma 
which follows the approach taken in |CJK1) and |C JK2j . 

Lemma 5.4. 



where C is a dimensional constant. 

Proof. 1. The proof is similar to that of Lemma 8.3 of |CJK2) . Note that u t is 
a caloric function in tt + (u) and Q~(u). Let us prove the lemma for u + , since 
parallel arguments apply to u~ . 

2. We divide ut into two parts. More precisely, let 



where both v\ and vi are caloric in f2 + (u), v\ has initial data zero and the 
boundary data a|-Du + | (|-Dt2 + | — |-D{t~|) on T(u), and V2 has the initial data 
u t (-, — 1/a) and the boundary data zero on T(u). 

3. As for vi, we need to use the absolute continuity of the caloric measure 
with respect to the harmonic measure, as well as the Lipschitz continuity of the 
free boundary, we proceed as in Lemma 8.3 of [CJK1| . Note that we have 

\Du+\ - |D?T| - 1 

in [-Bi/ 2 (e n ) U -Bi/ 2 (— e„)] x [—1/a, 1/a]: this follows from the assumption (|4.1[) , 
and Lemmas 13.31 and 13.41 Therefore we can proceed as in Lemma 8.3 of jCJKlj 
to obtain 



N < a\Du\ 2 < Ca in [B 1/2 (e n ) U B 1/2 (-e n )] x [-l/2o, l/2a], 



u t = vx + u 2 , 
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where w^ x ^ is the caloric measure for Q(u). 

vi(x,t) > a [ -|I>tr| ! W x '*> > -a\Du\ 2 (x,t). 

Jr(u)n{-i/a<s<t} 

4. As for V2, we conclude that it must be smaller than that of caloric 
function solved in the whole domain with the absolute value of its initial data. 
The advantage is that then we can use the heat kernel. Note that the initial data 
is given at t = —If a and has a compact support. The initial data is given by 
v t < ^v £ri , where v en (x, t) is comparable to the derivative of harmonic function 
in Lipschitz domain. 

Therefore the heat kernel representation is given as 



1 



Vn\ exp" 



-\x-y\ 2 /(t+l/a) v ( yj _l/ a )^. 



(f +1/o)t+ 

Since t € [0, 1/a], and fcexp _afe2 < Cexp~5 fe2 . we get the effect of O(a). 
Now we change the scale, and consider the function 



□ 



1 r 
v(x,t) := — u(rx + x ,— 1 + 1) (5.2) 
C r C 

Then this function is Lipschitz continuous, in space and time, away from 
the free boundary. The following lemma suggests that the cone of monotonicity 
improves away from the free boundary, as we look at smaller scales. The proof 
is parallel to that of Lemma 8.4 in \( S2|. 

Lemma 5.5. Let v given by (|5.2[) . Suppose that there exists constants 5 > 
and < A < B, fj, := B — A such that 

a(Dv, -e„) < S and A < < B 

-e„ ■ Dv 

in S 1 / 6 (— |e„) X (— Sf/j,, S/fj,) with — <r. Then there exist a unit vector v G M n 
and positive constants r ,&o < 1 depending only on A, B and n such that 

3 S 8 
a(Dv(x,t),f) < b S in B 1/8 (--e n ) x (-r -,r -). 

4 M t L 

Now we can proceed as in section 6 of jCJK2j to obtain further regularity, 
using Lemma [5.41 instead of the uniform upper bound on \Du\ up to the free 
boundary. 

Theorem 5.6. T(v) is C 1 in space in Qi/2- In particular, three exist constants 
Iq,Cq > depending only on L,n and M such that for a free boundary point 
(xo,to) £ r(w), T(v) n B 2 -i {xo, to) is a Lipschitz graph with Lipschitz constant 

less than if I > l . 
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Figure 4: Locally Lipschitz initial domain 
o Regularity in time 

Lastly, proceeding as in section 7-8 of |CJK2j yields the differentiability of 
T(v) in time. The main step in the argument is the following proposition: the 
statement and its proof is parallel to those of Theorem 7.2 in jCJK2j . 

Proposition 5.7. There exist constants Iq > and 1 < 7 < 2 depending only 
on L, 71, M such that for [xq, to) G F(w) fl Qi, if I > Iq then T(v) n B 2 -i (xo, to) 
is a Lipschitz graph with Lipschitz constant less than I 1 . 

Above proposition and the blow-up argument in section 8 of |CJK2j yields 
the desired result: 

Theorem 5.8. T(v) is differ entiable in time. Moreover 
C- 1 < \Dv\(x,t) <C inCl(u)nQ 1/2 , 

where C = C(M,n). 

6 General case: solutions with Locally Lipschitz 
Initial data 

In this section, we present how to extend the result of the main theorem to 
solutions with locally Lipschitz initial data. Our setting is as follows. Suppose 
f2o is a bounded region in Br(0). Suppose u is a solution of (ST2) with uq > — 1, 
1*0 = — 1 on -Br(O) and uq < Mq. Further suppose that Qq is locally Lipschitz: 
that is, for any xq G Tq, Tq n Bi(xq) is Lipschitz with a Lipschitz constant 
L < L n . 

Let the initial data uq solve Auq = in Bi(xq). Then we claim that the 
parallel statements as in Theorem 11.21 hold in £?2d„(a;o) x [t(xo, do)/ 2, t(xo, do)], 
where do is a constant depending on n and Mq. More precisely: 



34 



Theorem 6.1. Suppose u is a solution of (ST2) with initial data uq such that 
— 1 < uq < Mo. Further suppose that for xq <E Tq, Tq n B\(xq) is Lipschitz with 
a Lipschitz constant L < L n and Aug = in the positive and negative phases 
of uq in B\(xq). Then there exists a constant do > depending on n and Mo 
such that (a) and (b) of Theorem \1.2\ hold for u and d < do- 



The proof of the above theorem is parallel to that of Theorem 1 1.2 1 in section 
5, after proving the following lemma. 

Lemma 6.2. There exists a solution v of (ST2) with a star-shaped initial data 
such that the level sets of u and v are edo-close to each other in i?2d (* c o) U P 
to the time t(xo,do',u), where do > is sufficiently small. In particular, u and 
T(u) is (.-monotone in a cone of W x andWt in i?2d (xo) x [t(%Oi do)/2, t(xo, do)]. 

Even though our equation is nonlocal, the behavior of far-away region would 
not affect much the behavior of solution in the unit ball, if the solution behaves 
"reasonably" outside the unit ball. For example, in the star-shaped case, we 
know at least that the free boundary is almost locally Lipschitz at each time. 
In the locally Lipschitz case, we control the solution by putting an upper bound 
Mo on the initial data uq. We will argue that in a sufficiently small sub-region 
of Bi(xo) x [0, 1], the solution is mostly determined by the local initial data in 
Bi(xq). The perturbation method in the proof of Lemma 2.4 in jCJKl] will be 
adopted here. Denote Bi(xq) = B\. 

1. Construct a star-shaped region f2' C 5r(0) such that 

(a) n'nBj =Q n£i. 

(b) Q! is star-shaped with respect to every x <E K C ST for a sufficiently large 
ball K. 

Let Vq be the harmonic function in O' — K with boundary data 1 on dK, and 
on dD! . Next, let Vq be the harmonic function in -Br(O) — CI' with boundary 
data 1 on OBr(0), and on dft'. Let B2 be a concentric ball in B\ with the 
radius of e fc °, i.e., 

B-2 = B tkQ {xo) C B x {x ) = B x . 

Let ko be sufficiently large. Then by Lemma 12.61 a normalization of Vq by a 
suitable constant multiple yields that for any x € B2 

l- e <^4<l + e . (6.1) 
vq{x) 

Let v solve (ST2) with initial data w = Vq — . Then Theorem 11.21 applies for 
v since Vq is star-shaped with respect to K. 

For the proof of the claim, we will find a sufficiently small do such that v 
is edo-close to u in B2d ( x o) U P to the time t(xo,do). More precisely, we will 
construct a supersolution wi and a subsolution u>2 of (ST2) such that in some 
small ball Bh{xo), we have 

W2 < u < Wl 
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and the level sets of w\ and iu 2 are he close to the level sets of v. 

2. Let k\ and fc 2 be large constants which will be determined later. Define 
H± := (r Q (v)±e ko+kl e n )nB 2 . 

Let 

ct .— e 

and let t(do) : = i(£ojdo;i') = t(xo, do; u). First note that 

t(d ) > dl~ > e 7(fco+fc 1+ M/6_ 

Hence for w to be almost harmonic in a scale much larger than e fe °+ fcl ; we need 
V^o7> e fc °, i.e., 

7(fc + ki + fc 2 )/12 < ho. 
Observe that by the construction of H and do, 



Jt(d ) > radius(B 2 ) > distfiP , r ) > max distfx, r ) (6.2) 

x<£r t nB 2 ,o<t<t(d ) 

where the last inequality follows from Lemma 12.111 if we choose fc 2 > 2fci . If fc 2 
is sufficiently large, then one can prove from the last inequality of (|6.2[) and the 
bound on vt that 

l-e<-j — 7 rp = ~i -r-rr < 1 + e on i? 1 x 0,i(d ) ■ (6.3) 

l^oWl Fowl 

3. We do have an estimate, Lemma 12.111 on how far the boundaries move 
away for the local one-phase case. If we take the one-phase versions with initial 
data Uq and Uq , and compare with u, then we obtain that T(u) fl -B 2 stays in 

2-a 

the -neighborhood of T (u) D -B 2 up to the time i(do) = t(xo, do). In other 
words, the free boundary of u moves less than d^ 7 in B 2 up to the time t(do). 

Now we let S be the region between H + and H~. To construct a sub (or 
super) solution in S, we take the fixed boundary data (1 — e)vo(x) on H~ (or 
iJ + ), and (1 + e)vo(x) on H + (or ii" - ). To control the effect from the side 
dB 2 n S, we bend the free boundary T t (v) by d b J 7 on each side of dB 2 n 5, 
using the conformal mapping $ (or $). (See section 4 of for the definition of $ 
and $.) More precisely, we bend the free boundary of v downward (or upward) 
using the conformal map <t (or $), and solve the heat equation in there. Then 
similar arguments as in Lemmas 4.1 and 4.3 of [CK] yield that the solution is 
still (almost) a supersolution, and it stays close to the original solution. 
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(l+4r 2 )D t 




u>0 



4/ 3 



u<0 



